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Abstract.  The  problem  of  reconstructing  the  density  and  elastic 
properties  of  a slab  of  infinite  horizontal  extent  and  finite 
thickness  set  in  motion  by  a double  couple  excitation  is  considered. 
It  is  shorn  that  the  information  contained  in  the  amplitude  and 
frequency  response  associated  with  a single  horizontal  wave  number 
is  sufficient  to  insure  the  uniqueness  of  the  solution  of  this 
inverse  problem. 


Key  words : Inverse  problem  - Amplitude  and  frequency  response. 
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1.  Introduction 

The  present  paper  is  concerned  with  the  uniqueness  of  the  solution 
of  the  inverse  problem  for  the  internal  structure  of  the  earth. 

The  aim  is  to  find  what  data  are  required  in  order  to  find  both  the 
Lam6  parameters  and  the  density  uniquely.  These  data  should  be  such  as  to 
enable  us  to  perform  the  inversion  under  the  following  two  conditions: 

(i)  the  solution  should  be  obtained  ab  initio  and  should  not  be  a 
perturbation  of  an  existing  earth  model; 

(ii)  an  ad  hoc  assumption  & la  Adams -Williamson  should  not  be  used, 
but  rather  its  validity  should  be  checked. 

I have  been  able  to  make  some  progress  in  answering  the  question  of 
the  necessary  data  for  problems  simpler  than  that  for  the  earth  (Barcilon 
1976a,  1976b).  I believe  that  the  main  result  will  remain  valid  even  for 
the  more  complicated  geophysical  problem.  Namely,  if  the  source  is  known, 
then  the  amplitude  and  frequency  response  associated  with  a single  horizontal 
wave  number  is  sufficient  to  insure  the  uniqueness  of  the  inverse  problem. 

I shall  endeavor  to  present  an  overview  of  these  results  omitting 
most  of  the  underlying  calculations  which  either  have  appeared  or  will  appear 
elsewhere. 


3. 


1 


I 

iv  i 

I 

7 


2.  Slab  with  Double-couple  Excitation 

Let  us  consider  a perfectly  elastic  slab  of  infinite  horizontal 
extent  and  of  thickness  a_  (see  Fig.)*  We  shall  assume  that  the  density  p 

Insert  Figure. 

Caption:  Schematic  diagram  of  the  slab. 

and  the  Lam6  parameters  and  y are  solely  functions  of  the  depth,  i.e.  of 
jo  We  shall  also  assume  that  the  surface  x=0  is  stress-free,  i.e. 

=-«i  =%2  ' 0 at-'°' 

where  p..  stands  for  a component  of  the  stress  tensor.  The  slab  is  set  in 
motion  at  time  _t=0  by  means  of  a double-couple  applied  to  the  upper  sur- 
face £=<*.•  For  the  sake  of  presentation,  we  consider  the  simplest  case 
possible  and  write 


(2.1) 


at  x*a. 


% 3 

Exe, 3 ' T 6<£i>  a'«2>  Hw 

~~  ‘ T 6(£2)  HCt) 

In  the  above  formula,  6 stands  for  the  Dirac  delta  function  and  H for  the 
Heaviside  step-function;  a prime  denotes  differentiation  with  respect  to 
an  appropriate  variable. 

It  should  already  be  clear  that  the  model  we  are  considering  is 


(2.2) 


closely  related  to  the  ultimate  geophysical  problem.  The  slab  can  be 
looked  upon  as  the  mantle  of  a flat  earth  without  liquid  core  and  gravita 
tional  force.  To  pursue  the  analogy  with  the  geophysical  situation,  we 
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shall  refer  to  the  vertical  and  horizontal  displacements  u(a,C,t)  and 
va(a,£»t),  (ct=l,2),  at  the  upper  surface  as  the  seismograms.  The  in- 
verse problem  can  be  stated  thus:  given  these  seismograms,  can  we 

determine  a (x) , y (x)  and  p (x)  uniquely? 
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3.  Compressive  and  Torsional  Modes 

Just  as  for  the  earth,  the  response  to  an  arbitrary  excitation 

can  be  synthesized  by  means  of  a superposition  of  two  kinds  of  normal 

modes  of  oscillations  which  we  shall  refer  to  as  the  compressive  and 

torsional  modes.  In  fact,  in  the  spirit  of  Alterman  et  al  (1959),  let 

us  write  the  displacements  thus: 

00  00 

u(x,£,t)  » | e"---  du»  ||  YjCx.ujk)  H(£,k)  dk, 

• 00  00 

00  00 

VjOc.i.t)  = | e i—  doj  ||  {y3(x»u>;JO  jjr-  * j^r->dk, 

-00  00  ^ 

00  00 

“|  * d-  ||  W~  ■ ^"}d— 


In  the  above  representations 


H(C,k)  = exp(-i  C*k) 


is  the  function  appropriate  for  the  synthesis  of  the  horizontal  plan- 
form.  Thus  k»(k  ^3  is  a horizontal  wave  number  and  it  is  convenient 
to  define 

.2  .2  .2 

i.  ” - 1 + -2  ’ 

(u  stands  the  frequency  whereas  the  variables  y^  and  z ^ respectively 
are  associated  with  the  vertical  structure  of  the  displacements  occuring 
in  the  compressive  and  torsional  modes. 

In  addition  to  the  displacement  fields,  we  introduce  the  stress 
fields  as  follows: 


6. 


7. 


i'l 


0 


P+M 


1 

£ 

0 


with 


(3.7a) 


z_ 

= 0 

at  x=0 

-2 

•T 

(3.7b) 

i2 

COS  20=  0 

at  x=a. 

In  the  above  formulas,  (£, 6)  are  polar  coordinates  in  the  horizontal 
wave  number  space. 

We  shall  be  interested  in  the  above  equations  in  connection 
with  the  inverse  problem  for  which  we  shall  take  the  data  to  be 
y^(a,m;k),  y^CajWjk)  and  z^(a,u;k),  i.e.  the  temporal  and  spatial 
Fourier  transforms  of  the  seismograms.  This  statement  will  be  made 
even  more  precise  in  sequel. 
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4.  Properties  of  the  Seismograms 


Let  and  nj2^  (i=l,2,3,4)  be  two  fundamental  solutions  of 


(3.6a)  such  that 


_(!)  . & 

-l  lL 

nf2)  - 
-3  3j 


(j=l* 2,3,4)  for  x=0. 


(4.1) 


where  is  the  Kronecker  delta.  Clearly,  and  rij  ^ are  functions 

of  £,  (u  and  £.  The  solution  of  the  boundary  value  problem  (3.6)  can  be 
written  in  terms  of  these  two  solutions  as  follows: 

.iT  IVll?2'*  Ca,ou,  A) n (x,a>,£)+n^  (a,io,£)tj^  (x,uj,£)1 

y.  (x,co;£,9)  = -ii-  sin20  1 ± = — — — ± = ± 

“ 8jt  w “ W_ . (a,u;£) 

(j-1.2,3,4) 


(4.2) 


where 


W (£*“;£)  = 


n^1-*  (x,u>,£)  n^2^(x,u,£) 


jt^(x,u,£)  n^2-*  (x,u>,£) 


As  a result,  the  spatial  and  temporal  Fourier  transforms  of  the  seismo- 
grams have  the  following  forms: 


(4.3) 


k | 
. > 


iT  W12  CS,ju;  A) 

y (a,w;£,0)  = sin20  ^ f 
1 8ir  w 2!24la,a),£) 


y3(a,u;£,0) 


ir  W (a,u>;£) 

sin20  =7 rr-  . 

sA 


Looked  upon  as  functions  of  _m,  W^2(a,ai;£),  £23 and  W24(a,u;£) 
are  completely  determined  by  their  zeros.  Furthermore,  these  zeros  are 


(4.4a) 


(4.4b) 


mm 


natural  frequencies  of  compressive  modes  of  the  slab  in  three  vibrating 
configurations.  These  three  configurations  correspond  respectively  to 
the  following  boundary  conditions  at  x_=a; 

(i)  yj  = y2  3 ° 

(ii)  ll  = lz  3 0 

and 

(iii)  y2  = = 0 ; 

in  all  three  cases,  the  boundary  conditions  at  x_=0  are  = y^  = 0. 
Thus,  keeping  £ fixed  it  is  possible  (at  least  in  principle)  to  ex- 
tract two  "compressive"  spectra  over  and  above  the  compressive 
spectrum  corresponding  to  stress  free  conditions  at  x=0,  a. 

Similarly,  let  us  denote  by  ? (x,w;£)  a solution  of  (3.7a) 

—a.  — — 

such  that 


j 

(4.5a)  I 

j 

(4.5b) 


(4.5c) 


^(O.uj;*)  = 6la  (o=l,2). 


Then,  the  solution  of  the  boundary  value  problem  (3.7)  can  be  written 
thus 


zo(x,&i;*,e) 


cos20 


j^(a,wjS.) 


In  particular,  that  part  of  the  spatial  and  temporal  Fourier  transforms 
of  the  seismogram  associated  with  the  torsional  modes  has  the  form 


(a,o>;  £ , 0) 


it  ^1 

z — COS 20  ; :pr 

J>.. c„(a,uj;i) 
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Once  again,  the  zeros  of  C^(a,oj;£)  and  are  the  natural  fre- 

quencies of  the  torsional  modes  of  the  slab  in  two  vibrating  configura- 
tions corresponding  to  the  following  boundary  conditions  at  x=a 


(i)  zx  = 0 , 


z_2  = 0 . 


The  zeros  of  W^CajUjJi)  and  are  associated  with  the 
frequency  response  of  the  slab  whereas  those  of  ^12  (a»“;*)  > W23  (a,w;i.) 
and  c,  (a,a);£)  are  associated  with  the  amplitude  response. 
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5.  Solution  of  the  Inverse  Problem 


live  can  proceed  along  several  routes.  For  instance,  we  could 
work  with  the  so-called  influence  functions  for  a group  of  five  eigen- 
value problems  (Barcilon  1976a).  A second  approach,  particularly 
simple  for  discrete  versions  of  these  inverse  problems,  relies  on  con- 
tinued fractions.  This  method  has  been  pioneered  by  Krein  (1952)  and 
generalized  by  Barcilon  (1976b).  We  shall  refer  the  reader  to  these 
papers  for  a more  discursive  treatment  of  specific  problems.  Here, 
we  shall  confine  ourselves  to  the  derivation  of  the  formulas  essential 
for  carrying  out  this  procedure. 

The  first  step  consiscs  in  deriving  differential  equations  for 
the  six  distinct  Wronskians  IJL ^ (x,uj;£) . These  equations  are 


Jfc  = *2“l4-|%223- 

*i3  =1*14 

-14  = *12*t"“  *+4— 1 *13*  *24  * *34' 

*23  = *12  ' S*  *13  * 7 *24  - ‘2  *34’ 

“24  = -i-S  *14  * I***  * 4J!  ££  ^ *23’ 


*34 


-14  + £+2^  -23 


Actually,  from  (5.1a)  and  (5. If)  together  with  (4.1)  we  deduce 
that  W34  and  are  linearly  related,  viz. 

W12(x,(o;£)  + £2W34  (x,oj;  £)  = 0. 


12. 


The  Wronskians  satisfy  also  the  following  identify 

— 12  «34  * —13  ^42  * -14  *23  ' ° • 

or  in  view  of  (5.2) 

J^2(x)  + ^2W13(x)W24(x)  - *2W14(x)W23(x)  = 0. 

By  eliminating  also  from  (5.1)  and  incorporating  the  equations  (3.7a) 
for  we  obtain  the  following  grand  system  of  first  order  equations: 


0 0 V 


-2k  2 413L,2 


A+2v. 


1 

2.+2£ 


1 

X_+2^ 


0 “13 


0 


0 —23 


2 2 ^+u  .2  n 

-a.  £ e.+4Hjfe2u  0 


0 0 w24 


u U £ -1 
0 -o;2£+viJ12  0 £2 


If  the  slab  is  discretized  in  N layers,  then  W.^(a,w)  and  ^(a.w)  are 

2 “ " 

known  polynomials  of  degree  N in  jo  . Indeed,  W12(a,oj),  W^fa,^)  , ' 

£^(a_,<j)  and  ^2(a,w)  are  extracted  from  the  seismograms,  whereas  W^(a,io) 
and  Wj,4(a,oj)  are  deduced  from  (5.3).  Furthermore  the  W' s and  ?'s  are  also 
known  at  x=0,  viz 
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